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Abstract. We present a detailed formalism for the description of collective electron excitations in fullerenes 
in the process of the electron inelastic scattering. Considering the system as a spherical shell of a finite 
width, we show that the differential cross section is defined by three plasmon excitations, namely two 
coupled modes of the surface plasmon and the volume plasmon. The interplay of the three plasmons 
appears due to the electron diffraction of the fuUerene shell. Plasmon modes of different angular momenta 
provide dominating contributions to the differential cross section depending on the transferred momentum. 
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1 Introduction 

In the present paper, we give a detailed theoretical ex- 
planation of the formation of various plasmon excitations 
in fullerenes caused by collision with fast electrons. We 
demonstrate that a non- uniform electric field of a charged 
projectile causes variation of the electron density on the 
inner and the outer surfaces of the molecule as well as 
the volume density variation inside the fuUerene shell. The 
variation of electron density leads to the formation of three 
plasmon excitations, namely two surface plasmons and the 
volume plasmon. 

Plasmon excitations represent the collective oscillation 
of electrons of an atomic system against the positively 
charged ions. This collective electronic motion appears 
when the system is exposed to an external field (the elec- 
tromagnetic field or the electric field of a charged projec- 
tile). Collective excitations manifest themselves in the for- 
mation of giant resonances in the excitation spectrum of 
atomic systems. The existence of plasmon resonances is a 
general phenomenon occurring in various atomic systems, 
while a position of the resonance depends strongly on the 
type of the system. For instance, collective excitations in 
many-electron atoms have typical resonance frequencies of 
about 100 eV [1]. Plasmon excitations in various atomic 
clusters have much lower frequencies, namely of the order 
of several eV for metal clusters [2,3] and of several tens 
eV for fullerenes [4,5]. 

Delocalized electrons of atomic clusters may form two 
different types of collective excitations, namely the sur- 
face and the volume plasmons [6,7]. The dipole surface 
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plasmon is responsible for the formation of the giant res- 
onance in photoabsorption spectra of metal clusters [8, 
9] and fullerenes [4], and play also an important role in 
the process of inelastic scattering of electrons [10, 11]. The 
excitation of metal clusters revealed the existence of the 
volume plasmon which has a higher resonance frequency 
and is essential for the formation of the electron impact 
ionization cross section [12]. 

Existence of the giant resonance in the excitation spec- 
tra of fullerenes at about 20 eV was predicted theoretically 
[13] and then observed experimentally for the gas phase 
Ceo compounds studying the processes of photoioniza- 
tion [4] and inelastic scattering of electrons [5] . Formation 
of the dipole plasmon resonance in the photoionization 
cross section of Ceo was studied theoretically within var- 
ious models and approaches [14,15]. Recent experiments 
on photoionization of neutral [16] and charged [17] Cgo 
molecules revealed the existence of the second collective 
resonance at about 40 eV which later was associated [18] 
with the second surface plasmon. 

Theoretical investigations of the scattering of fast elec- 
trons on fullerenes [19,20] revealed the existence of the 
diffraction phenomena arising in the scattering processes. 
The first experimental observation of the electron diffrac- 
tion on Ceo was reported in Ref. [21]. It was shown that 
plasmon modes of different angular momenta provide dom- 
inating contributions to the differential cross section at 
different electron scattering angles. In the cited papers, 
the fuUerene was modeled as an infinitely thin sphere and 
collective electron excitations were represented by the sin- 
gle surface plasmon. 

In Ref. [22] a hydrodynamic model was applied to de- 
scribe the plasmon excitations formation mechanism. It 
was shown that no volume plasmons can be excited in 
the dipole-photon limit. Therefore, the volume plasmon 
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can manifest itself only when the system interacts with a 
non-uniform external field, e.g. in collisions with charged 
particles. 

The paper is organized as follows. Section 2 is dedi- 
cated to the detailed description of the theoretical frame- 
work. We introduce the model of a fullerene in Sec. 2.1. 
In Sec. 2.2.1, we introduce general expressions based on 
classical electrodynamics and hydrodynamics describing 
collective electron excitations in a many-electron system. 
In Sec. 2.2.2, we derive a general expression for the elec- 
tron density variation in an arbitrary spherically symmet- 
ric system exposed to an external field. In Sec. 2.2.3, we 
apply this formalism to a fullerene considered within the 
introduced model. In Sec. 2.3, we derive general expres- 
sions for the inelastic scattering cross section. In Sec. 2.4, 
we present the expressions for the volume and surface plas- 
mon excitations obtained within the plasmon resonance 
approximation. We discuss also the diffraction phenomena 
which manifest themselves in the contribution of terms of 
different multipolarity. In Sec. 3, we consider two limit- 
ing cases of the general equations presented in Sec. 2. We 
consider the "metal cluster" limit when the system is rep- 
resented by a charged full sphere, thus the volume and the 
single surface plasmons appear in the system. We consider 
also the "infinitely thin fullerene" limit when collective 
excitations are described only by the single surface plas- 
mon. The obtained expressions coincide with the results 
presented earlier in Refs. [12,19,20]. In Sec. 4, we draw 
the conclusions from this work. In Appendix, we describe 
the transformation of general expressions obtained in Sec. 
2 in the case of a uniform external field. We show that 
no volume plasmon can appear in the system exposed to 
a uniform field and the photoionization cross section is 
defined then by the two surface plasmons. 

The atomic system of units, nie = e = h = 1, is used 
throughout the paper. 

2 Theoretical framework 
2.1 Model of a fullerene 

In this paper, wc consider a fullerene Cn as a spherically 
symmetric system where the charge is distributed homo- 
geneously between two concentric spheres [23-25]. The 
width of the fullerene is AR = i?2 — i?i where Ri , i?2 are 
the inner and the outer radii of the molecule, respectively. 
The equilibrium electron density distribution po{r) is ex- 
pressed via the number N of delocalized electrons (four 
2s^2p^ valence electrons from each carbon atom) and the 
fullerene volume V: 

_ ( N/V ioiRi<r<R2 

\ if otherwise . ^ > 

The volume of the fullerene shell can be expressed as 

V=f{Rl-Rl) = '^Rl{l-e), (2) 

where ^ = Ri /i?2 < 1 is the ratio of the inner to the outer 
radii. 



This model is applicable for any spherically symmetric 
system with an arbitrary value of the ratio ^. Supposing 
^ = 0, one obtains a model of a metal cluster, while the 
case ^ = 1 represents a fullerene modeled as an infinitely 
thin sphere. 

When the fullerene is exposed to an external field, the 

electron density is modified, p = po+Sp. Below, we present 
a detailed formalism of the electron density variation, 6p. 

2.2 Variation of the electron density 

2.2.1 General equations 

Let us describe briefly a simplified version of the approach 
introduced in Ref. [22]. The presented formalism describes 
the behavior of collective electronic excitations and is based 
on classical electrodynamics and hydrodynamics. 

Let /9o(r) denote the stationary distribution of the elec- 
tron charge in the point r. The variation of the electron 
density, Sp{r,t), depends on the position r and time t. 
Therefore, the total electron density is introduced as: 

p{r,t)= po{r)+Sp{r,t) . (3) 

Following [22], let us describe the collective motion 
of the electron density using the Euler equation and the 
equation of continuity. 

The Euler equation couples the acceleration, dv(r, t)/dt, 
of the electron density with the total electric field E acting 
on the system at the point {r,t): 

^=E(M). (4) 

The electric field E includes both the external field acting 
on the system and the polarization contribution due to 
the variation of electron density 5p{r,t): 

where (f>(r, t) is the scalar potential of the external field. 

Introducing (5) in (4) and evaluating the full time 
derivative of the vector v one obtains 

2xM+(v(r,,).v)v(,.,t) = 

-(v«M))-v/dr'^£<r:i^^. (6) 

The potential of the external field is assumed to satisfy 
the wave equation and has the monochromatic dependence 
on t: 

0(r,t)=e-*0(r) , (7) 
where 4>{r) satisfies the equation 

Zi0(r) + 0(r) = (8) 

with q is the wave vector. 
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The motion of electron density in the system obeys the 
equation of continuity, which reads 



^ + V.(p(M)v(M))=0. 



(9) 



Eqs. (6) and (9), being solved simultaneously, deter- 
mine the variation of electron density 5p{r, t) as well as 
its velocity v(r, t). 

Let us estimate the relative value of the first and the 
second terms on the left-hand side of Eq. (6). To make 
an estimation, let us consider the one-dimensional case, 
and let E oc cxp(iti;<:). Thus, Eq. (6) can be written in the 
simplified form: 



ox 



Here dv/dx ^ v/R where R is the size of the system. 
Supposing that the term i) dominates, one finds that v ~ 
(l/cij)i?o cxp(icjt). Hence, from the condition i) ^ v'^/R 
one can obtain the criterion E <C uj^R. Below, we assume 
that this condition is fulfilled and neglect the second term 
on the left-hand side of Eq. (6), which means physically 
that the external field causes only a small spatial inho- 
mogeneity in the electron density distribution within the 
system. 

Let us seek the solutions of Eqs. (6) and (9) in the 
following form: 



6p{r, t) = 5p{t) e''^* 
v(r,t) = v(r)e'"* . 



(11) 



Substituting these expressions into Eqs. (6) and (9) 
and performing some transformations with the simultane- 
ous use of Eq. (8) and A\t — r'|~^ = — 47r(5(r — r'), one 
derives a set of the following linear equations: 



v(r) 



1 



V(t>{r) + V 



/ 



5p{v') 



dr' 



(12) 



and 



{uj^ - 47rpo(r))5p(r) + Vpo(r) • V j dr' 

= g2po(r)^(r)- Vpo(r)-V</.(r) . 



(13) 



2.2.2 Arbitrary spherically symmetric system 



Eqs. (12) and (13) describe the dynamics of electron den- 
sity under the action of an external field. In the case of the 
spherically symmetric density distribution, po(r) — Po{r), 
one can exclude angular variables from Eq. (12) and (13). 

Let us expand functions ^(r), 5p{v) and |r — r'|~^ into 
spherical harmonics: 



'(/.(r) = ^(/),(^)i1™(r) 

Im 

5p{v) = Y,hi{r)Yira{y) 

Im 

Tr37^ = E^^'(^'^')^™(r)F4(r'), 



(14) 



where the function bi{r,r') is defined as follows: 



e{r'-r)+^^0{r-r') , 



0{x) is the Heaviside step function defined as 



x<0 




(15) 



(16) 



and the notation ilj = y/2l + 1 is introduced. 

Using the last two expressions in Eq. (14), one gets 



(10) / ^dr' = Y: ^YUr) J r'%iry)Sp,ir')dr' . 



(17) 

To expand the terms containing the operator V, the 
following general expression is used (see Ref. [26]): 



'(/(r)F;„(r)) = 



_ d/(r) v(-i) 



(r) 



+ V^(«TT)^Y«(r), (18) 

where Y^^^^y) and Y[^^(r) arc the longitudinal and the 
transverse vector spherical harmonics, respectively. The 
definition of these functions can be found in Ref. [26]. 
Hence 

oo 



dbiir, r') 
dr 



(r) + 7^M^YL)(r) 



(19) 



and 



V</>(r) = ^ 



+ 



^AaTT)^Y«(^) 



Accounting for the identities (see Ref. [26]) 

YL;/'(r)-Y«,(r)=0 
YL''(r)-Y,tV(r) = 11„(r)W(r), 
Eq. (13) is transformed to the following expression: 



(20) 

(21) 
(22) 



E^™(^ 



-I- 47r 



-A'KpQ{r)j 6pi{r) 

OO 

Po(0 



gi{r, r') Spi{r')dr 



Im 



Y,yim{r)[q^Po{r)Ur)-p'o{r)m] , (23) 



Im 
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where the function gi{r,r') is related to bi{r,r') defined in 
(15): 



gi{r,r') = r' 



,2 dbi{r,r') 



dr 



(24) 




Multiplying both sides of Eq. (23) by Y{!^{y) and carry- 
ing out the integration over spherical angles of the vector 
r, one obtains a general equation for the variation of elec- 
tron density in an arbitrary spherically symmetric system: 



Fig. 1. Left panel: Representation of a fuUerene as a spherical 
shell of a width R2 — Ri. Variation of the surface charge den- 
sities, cr'^'^', and the volume charge density, 6q, is also shown. 
Right panel: Representation of the symmetric (a) and the an- 
tisymmetric (b) modes of the surface plasmon. 



00 

/ \ p' (r) f 

{^u"^ — 4:TTpo{r)jSpi{r) + 4tt J gi{r,r')5pi{r')dr' where the following notations are introduced: 



q^Po{r)Mr)-Po{rmr) 



(25) 



where <j)i (r) is the scalar potential and co is the frequency 
of the external field. 



2.2.3 Fullerene as spherical shell of a finite width 

In the case of the fullcrcno model described in Sec. 2.1, 
the equilibrium electron density distribution, po{r), is con- 
stant within the interval Ri < r < R2 and equals to zero 
if otherwise: 

Poir)=PoO{r-Ri)e{R2-r) , (26) 

where po is defined by Eq. (1). The derivative of the func- 
tion po{r) is given by 

p'oir) = po [6{r - R,) - 6{r - R2)) , (27) 

where S{x) is the delta-function. 

Then, the solution of Eq. (25) for the fullerene is sought 
in the following form: 

Spiir) = Sgi{r)e{r ~ i?i)©(i?2 - r) 

+ ap^(5(r - i?i) + a^Sir - R2) , (28) 

where SQi{r) describes the volume density variation aris- 
ing inside the fullerene shell, and crp'^-* are variations of 
the surface charge densities at the inner and the outer sur- 
faces of the shell, respectively (see the left panel of Fig. 
1). The volume density variation causes the formation of 
the volume plasmon, while the variations of the surface 
densities correspond to two surface plasmon modes. 

Using (26), (27) and (28) in Eq. (25) and carrying out 
algebraic transformations, one derives: 

{w - l)Sgi{r)e{r - Ri)0{R2 - r) 



(2) , r(2) 



t2 "l S 



n, 



]d{r-R2) 



e{r - Ri)0{R2 - r) 



— [<l>'i{Ri)6{r - i?i) - <P[iR2)S{r - R2] 



(29) 



w = 



(30) 



R-2 

^i-i f ^eijx ) 



da; 



7-(2) _ 
1 — 



1 1 



R2 



(31) 



Ri 



Parameter (jjp IS the volume plasmon frequency associated 
with the density po- Neglecting the dispersion, the volume 
plasmon frequency has a constant value and is defined as: 



2 . 3iV 

^p=47rpo=^3(i_^3) 



(32) 



Matching the terms of different types on the right- and 
the left-hand side of Eq. (29), one obtains three equations: 
one for the volume plasmon and the other two - for the 
surface plasmons. 

The solution corresponding to the volume density vari- 
ation reads as: 



Sgi{r) 



Mr) 



(33) 



IS 



w — 1 An 

The total density variation due to the surface plasmons 
ai{r) = ap^5(r - i?i) + ^(r - R2) , (34) 



where the quantities aj^^ and crp^ satisfy the following 
system of coupled equations: 



w ■ 



l + l 



.(1) 



1_(2) 



n. 



2? ^- 



(2) 



F2 



where 



Fi,2 = T 



mi 



An 



-I 



(1>2) 



(35) 



(36) 
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Using the expression (33) for the volume density vari- 
ation, the functions ij^^ and can be rewritten in the 
following form: 



R2 



n 



(1) 



w 



I f R[-' Mr) 



r(2) 



i + i 



47r 



da; 



w 



R2 

I 



J+2 



Mr) 



(37) 



1 nf J 47r 

Ri 



dx . 



The determinant of the system (35) is 

A= {w- Wil)(w - W2l) , 



(38) 



where wu and W21 are the roots of the secular equation 
A = 0: 



Wll 



W21 



1 



2 V 21 + 1 

1 



21 + 1 



l + 4i(Z + l)^2(+i^ 
1 + 4l{l + l)^2J+i ^ 



(39) 



Solutions of the system (35) are given by the following 
expression: 



/I) - 



1 



FiW-j^{Fi+e-'F2) 



A 



F2W-'-^{F2+e+'F,) 



T2 
A 



(40) 



Variation of the surface charge densities, crp'^\ results 
in the formation of two coupled modes of surface plasmon 
oscillations [23 25]. Frequencies of the symmetric, wi;, and 
the antisymmetric, u^i-, surface plasmons of multipolarity 
I are given by the expression [24,25]: 



^jl=WjiU)l (j = l,2). 



(41) 



In the symmetric mode the charge densities on the two sur- 
faces oscillate in phase, while in the antisymmetric mode 
they are out of phase (see the right panel of Fig. 1). 

Using Eqs. (33) and (40) in (28), one obtains the ex- 
pression which defines the multipole variation of electron 
density in a spherically symmetric fuUerene of a finite 
width under the action of the multipole component (r) 
of the external field: 



6pi{r) 



q 



Mr) 



w ■ 



1 An 



e{r-Rt)G{R2-r) 



+ ^5{r-R,) + ^5{v 



2.3 Inelastic scattering of an electron 



Ri 



(42) 



In the process of inelastic scattering the projectile electron 
undergoes the transition from the initial electron state 
(£i,Pi) to the final state (e2,P2) which is accompanied 



J 




I 



Fig. 2. Diagrammatic representation of the inelastic scattering 
process. The projectile electron goes from the initial electron 
state (ei,pi) to the final state (e2,P2) while the fuUerene is 
ionized from the initial state i (sj) to the final state / (e/). 



by the ionization (or, excitation) of a fuUerene from the 
initial state i with the energy to the final state / with 
£/. Diagrammatic representation of the process is given in 
Fig. 2. 

The matrix clement, M, which defines the amplitude 
of the inelastic scattering is given by 



M 



/,2 



1 



1,1 



= E/ ^^^*(r)V'/(K})| 



X Vi({ra})Vi+Hr){drJdr , (43) 

where {va] = i"i . • . tat are the position vectors of the de- 
localizcd electrons in the fullercne, r is the position vector 
of the projectile, ^'\^\y) and -02 \'^) stand for the initial- 
and the final state wave functions of the projectile, respec- 
tively. Superscripts {+) and (— ) indicate that asymptotic 
behavior of the wave functions is 'plane wave -|- outgo- 
ing spherical wave' and 'plane wave -|- incoming wave', 
respectively. 

The matrix element can be written as follows: 



M 



I 



An dq 



(2 I 



-iqr 



1) / 



(44) 



where q = pi — p2 is the transferred momentum. 

In the present paper, we consider the collision process 
of Ceo molecules with fast electrons. Since the collision 
velocity is larger than the characteristic velocities of de- 
localized electrons in the fuUerene, the first Born approx- 
imation is applicable [19]. Within this approximation the 
initial and the final states of the incident electron can be 
described by plane waves: 



^ipi-r 



4~ 



(r) 



(45) 



Within the framework of the plane-wave first Born ap- 
proximation the amplitude of the process reduces to 




(46) 



q=pi-P2 



The magnitude of is related to pi^2 and the scattering 
angle 9 = prP2 via: 

q'^=pl+pl- 2piP2 cos e « pl9^ . (47) 
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The final approximate equality is valid when pi w p2 and 
the scattering angle is small, ^ <C 1 rad. 

Performing the multipole expansion of the exponential 
factors in (46) (see, e.g., [26]), one obtains: 



M = 47r^i'y,;,(q)(/ 



where we introduced the following notation: 
Mr) =4.'^ 



(48) 



(49) 



and ji is the spherical Bessel function. 

Let us consider a general expression for the cross sec- 
tion of the scattering process: 



d. = ^^(.,.-.)EEl^l^7|y^dp, 

pol^ pol^ 



(50) 



where w = £i — £2 is the energy transfer, ojfi = Sf — Si 
and Lo = ujfi due to the energy conservation law. The sign 
Spoij, denotes the summation over the projection of the 

final state / orbital momentum, whereas X)poi denotes 
the averaging over the projections of the initial state or- 
bital momentum, dpf is the density of the fuUerene final 
states. 

Substituting the scattering amplitude (48) into Eq. 
(50), one derives the triply differential cross section: 



d^a 



d£2di7p 



1P2 



E 

Im 



pol 



xd{cOfi - co)dpf , (51) 



where 



(52) 



is the multipolar potential of the projectile electron, di^p^ 
denotes the differentiation over the solid angle of the scat- 
tered electron and sign J dpf means the summation over 
the final states (which includes the summation over the 
discrete spectrum and the integration over the continuous 
spectrum). 



2.4 Plasmon resonance approximation 

Let us consider the inelastic scattering cross section within 
the plasmon resonance approximation [19,20]. It relies on 
the assumption that collective plasmon excitations give 
the main contribution to the cross section in the vicinity of 
the giant resonance. Hence, one can neglect single-particle 
excitations when calculating the matrix element in Eqs. 
(50) and (51). 

According to the Kubo linear response theory [19,27], 
the integral on the right-hand side of Eq. (51) can be re- 
lated to the density variation 6pi (42), so one can perform 



the following substitution: 



5{u>fi- u>)dpf 



pol 



{r)Spi{uj,q;r)dx . (53) 



Here 6pi{oj,q;r) is a partial density variation due to the 
exposure of the system to the multipolar potential Vim (r) . 
In a general case, this variation depends on the frequency 
CO, transferred momentum q and the position vector r as 
well. 

Using (53) in (51), the triply differential cross section 
acquires the form: 



dV 



d£2di7p2 



J_P2 
TT^ Pi 



Elm[/;(u;,g)] , (54) 



where 



Ii{u;,q)=nfJ V;^ir)6pi{u;,q;r)dr 



(55) 



and the notation 77; = ^21 + 1 is used. 

Using the multipole variation of the electron density 
Spi{uj,q;r) defined by Eq. (42) as well as the notations 
(49) and (52), one can write 



Iliu, q) = q) + I\''>{i^, q) + ll''>{w, q) , (56) 

where 



00 

itH^'l) = j r'^ji{Qr)Sgiiuj,q;r)dr 



00 

ll''\u;,q) = ^nfj r^ji{qr)al^\co,q;r)dr , 



(57) 



and j = 1,2. The term 6Qi{u),q;r) is the electron den- 
sity variation associated with the volume plasmon (see 

Eq. (33)), a^^^ and ap' describe the density variation at 
the inner and the outer fullerene surfaces (see Eq. (40)). 

Performing some transformations, we come to the for- 
mula for the differential inelastic scattering cross section 
with no damping of plasmon oscillations: 



dV 



2^P2 

d£2di7p2 TT Pi 

r,.,2 

X imy~^ 



"21 



(58) 



where ojp is the volume plasmon frequency defined by Eq. 

(32) , ujii and UJ21 are the frequencies of the symmetric and 
antisymmetric surface plasmons of multipolarity I defined 
by Eq. (41). Functions Vi{q), Su{q) and S2i{q) are defined 
as follows: 
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Carrying out the imaginary parts produces: 



Vi{q) = nf zl (^jf (^2) + if (^2) - e (jfizi) + j^z,)) 

+ Sz2(^jl{z2)j'i{z2) - ^'^jl{zi)j'i{zi)^ 

-l[l + l)(jf{z2)-UKzi)) 

9772 / >^ 

- ((1 - «2)if (-^2) + ^(1 - ai)if (^i)) 



+ 



3i{zi)3i{z2) 



1 _ ^2;+i 



(59) 



and 



5u(g) = 2il; 

S2i{q) = 2ni 



02 - wii)ji{z2) + 

(1 - ^2i+l)(a2 - Wu){W2l - wu) 

(y{w2l - a2)jl{Z2) - ^'■~^'^W2ljl{zi)^ 
(1 - $^'+^)(W2i - a2){W2l - Wii) 



(60) 



where zi^2 = qRi,2, wu and W21 are defined in Eq. (39) 
and the following notations are used: 



ai 



(12 



l + l 

n? 



(1 - ^2'+!) . 



(61) 



The functions Vi{q)^ Sii{q) and S2i{q) are the diffrac- 
tion factors depending on the transferred momentum q. 
They determine the relative significance of the multipole 
plasmon modes in various ranges of electron scattering an- 
gles. Dominating contribution of different multipole modes 
results in the significant angular dependence of the differ- 
ential electron energy loss spectrum [28]. This phenomenon 
arises due to the electron diffraction of the fullerene shell 
[21]. Considering a fullerene as an infinitely thin sphere, 
the dominating contribution of different multipole modes 
to the surface plasmon was shown in Refs. [19,20]. 

Since plasmons decay from the collective excitation 
mode to the incoherent sum of single-electron excitations, 
it is essential to account for the damping of plasmon oscil- 
lations. Therefore, one should introduce the finite widths, 

rj-^^ and Fji ^ (j = 1, 2) of the volume and surface plasmon 
resonances, respectively, and make the following substitu- 
tions in the right-hand side of Eq. (58): 



1 



1 



a;2 - a;2 



jl 



a;2-a;2 +iwrif) 



up- - W2 



;2 - a;2 + -^^A^) 



(62) 



Im 



(s)2 
3I 



Im 



(v) 



w2 _ w2 + iwr/"' (a;2-w2)2^^2/^M2 



.(63) 



Thereby, the final formula for the differential inelastic 

scattering cross section with the account for three plas- 
mons and with the damping included is the following: 



d^a 



JS^W d3^(.i) d3^(.2) 



de2dJ7p2 
where 



d£2di7p2 d£2dl2p 



d£2di7p 



(64) 



d3^(-) 


_ 2i?2P2 








Vl{q) 


de2di7p2 


T^q^vi 


(w2 


--If 




dV^'i) 


2R2P2 


-E 

I 






Suiq) 


d£2di7p2 


nq'^pi 


(w2 








2R2V2 


-E 




,,2 p(s) 
-21 ^ 21 


S2i{q) 


d£2di7p2 


nq^pi 


(u;2 


--lit 


+ c.2rif 



(65) 



and the functions Vi{q), S\i{q) and S2i{q) are defined above. 

Let us define the angular momentum range which is 
considered with the introduced model. In Ref. [19], it was 
shown that excitations with large angular momenta I have 
a single-particle nature rather than a collective character. 
It follows from the fact that with increasing I the wave- 
length of the plasmon mode becomes smaller than the 
characteristic wavelength of the delocalized electrons in 
the fullerene [19]. In the case of the Cgo fullerene, only 
terms with I < 3 should be included to the sum over I in 
Eq. (65), while multipole excitations with Z > 3 are formed 
by single-electron transitions. 

The introduced model is applicable within the long 
wavelength limit, when the characteristic scattering length, 
1/q, is large. Under the condition of the small transferred 
momentum q, the volume plasmon is characterized by the 
constant frequency ujp which does not depend on the trans- 
ferred momentum [29]. In this paper, we do not consider 
the dependence of the plasmon widths on the transferred 
momentum which was studied in Ref. [12]. The widths are 
treated as external parameters which are not calculated 
within the present model. 



3 Limiting cases of the general formulae 

General expressions (64) and (65) for the differential in- 
elastic scattering cross section are applicable for any spher- 
ically symmetric system with an arbitrary value of the ra- 
tio ^. In this section, we consider a transformation of the 
general expressions in the two limiting cases: a metal clus- 
ter (^ = 0) and a fullerene modeled by an infinitely thin 
sphere (^ = 1). 
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3.1 The 'metallic cluster limit' 

Metal cluster can be treated as a system with a uniform 
electron density distribution over the sphere of a radius R 
[30]: 



N 



(66) 



where V = A-kR^ /2> is the cluster volume and N is the 
number of delocalized electrons in a cluster. In this limit, 
one obtains ^ 0, i.e. i?i — > and i?2 = R- 

Electron density variation on the cluster surface leads 
to the formation of the surface plasmon, while the volume 
plasmon arises due to the density variation inside the sys- 
tem. In this case, the antisymmetric surface plasmon mode 
does not contribute to the cross section: 



hm ■ = 

5->-o de2dJ/] 



(67) 



P2 



and the general expressions (64) and (65) for the differen- 
tial cross section transform into the following: 



d^a 



d£2dl2p2 d£2di7p2 d£2dJ7p2 



where 



d3(7(^) _ 4i?p2 'f- nftofrl'^jfiqR) 



d£2d/2p2 7rg4pi ^ ^^2 _ ^ ^2 JJ' 



,(s)2 



d3cr(^) 2R^P2 V- 



d£2d^2p2 Wpi ^ (^2 - + w2r; 



(v)2 



(68) 



(69) 



(70) 



and 



F{^) = jf{z) - ji-i{z)ji+i{z) - -ji{z)ji+,{z) . (71) 

z 

Within this limit, the general expressions for the frequen- 
cies of the volume (32) and the surface (39) plasmons 
transform into the following: uip = \J^N/R^, 

UJl = ^/l/{2l+l)LUp. 

The differential cross section (68) of the inelastic scat- 
tering on a metal cluster defined by Eqs. (69) and (70) 
coincides with expressions presented earlier in Ref. [12]. 



Therefore, the general expressions (64) and (65) reduces 
to: 



d£2dl2p2 irq^pi Y (^2 - wf^ % uj^r^ 
where 



' l{l + l)N 
(2Z+l)i?3 



(73) 



(74) 



is the surface plasmon frequency, and /) = is its 

width. Eq. (73) coincides with the expression presented 
earlier in Refs. [19,20]. 



4 Conclusion 

In this paper, we have presented a detailed formalism for 
plasmon excitations in fuUerenes caused by the collision 
with fast electrons. We have demonstrated that the energy 
loss spectrum is formed by three contributions, namely 
two modes of the surface plasmon and the volume plas- 
mon. 

Within the presented model, a fullerene was considered 
as a spherical shell of a finite width where the negative 
charge was distributed homogeneously over the shell. We 
showed that exposure of the system to an external field 
of a charged projectile causes the formation of the surface 
electron density variation on the surfaces of the shell and 
the volume density variation inside the shell. These vari- 
ations lead to the formation of three plasmon excitations. 

The introduced model can be applied to a spherically 
symmetric system with an arbitrary width of the shell. 
We have considered two limiting cases of the model when 
the system is treated as a charged full sphere ('metal clus- 
ter' limit) and as an infinitely thin sphere. The obtained 
expressions coincide with the results obtained earlier in 
Refs. [12,19,20]. 

Numerical results of the calculation based on the pre- 
sented formalism as well as the comparison with the recent 
experimental data will be presented in another paper of 
this issue of the journal [28] . 

A.V.V. is grateful to Deutscher Akademischer Aus- 
tauschdienst (DAAD) for the financial support. 



3.2 The 'infinitely thin fullerene' limit 

The model of a fullerene as an infinitely thin sphere was 
applied in a number of papers studying the photoioniza- 
tion [14] and the electron scattering [19,20] processes. In 
this limit, one obtains Ri ^ R2 = R and C 1- 

In the case of a sphere, the volume plasmon and the 
antisymmetric surface plasmon mode do not contribute to 
the cross section: 

5^1 d£2dl?p2 d£2dr2p2 



A Interaction with a uniform external field 

In this Appendix, we consider the transformation of gen- 
eral expressions for the electron density variation obtained 
in Sec. 2.2.1 and 2.2.3 in the case of the uniform external 
field. This case describes the interaction with an electro- 
magnetic field. 

We assume that the wave length of electromagnetic 
radiation is much larger than the typical size of the sys- 
tem, i.e. the condition ujR ^ 1 is fulfilled. This condition 
implies the validity of the dipole approximation. 
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In this limit, q = and Eq. (25) for the multipole 
variation of the electron density in a spherically symmetric 
system turns into the following one: 

(iv^ - 4npo{r))5pi{r) + j gi{r,r')5pi{r')dr' 



Therefore, one obtains the following expression for the 
total variation of electron density due to the surface plas- 
mons: 



A 4tt 



w ■ 



l + l 

m 



(i-e 



5{r - R2) 



(75) 



- yw-jp{l-e-')]S{r-Rr) 



(82) 



Carrying out the transformations described in Sec. 2.2.3 
and taking into account that in the dipole approximation 
the field intensity does not depend on the spatial coor- 
dinate 

0;(i?i) = 0j(i?2) = 0; , (76) 

one derives the following equation: 
{w - l)5Qi{r)0{r - i?i)e(i?2 - r) 



For the further discussion let us introduce the multi- 
pole moment, Qi, induced by the external field 0;: 



(1) 



(2) 



Ml 



(77) 



Matching the terms of different types on the right- and 
the left-hand side of Eq. (77), one obtains the equations 

for the volume and the surface electron density variations. 
The equation for the volume density variation reads: 

{w - l)5Qi{r)e{r - Ri)0{R2 - r) = . (78) 

Therefore, there is no volume plasmon arising in the sys- 
tem due to interaction with the uniform external field. 

The total electron density variation is described only 
by the surface contributions: 



Spi{r) = ap'(5(r - Ri) + cjl^'5{r - R2) 



(2), 



(79) 



which can be obtained solving the following system of 
equations: 



w ■ 



'+2.^(1) 



>S "l 



I 



w - 



nf J^' 



(80) 



(2) 



An 



The solutions of Eq. (80) are 



i + i 



where A — {w — w\i){w — W21) is the determinant of the 
system (80). 



Qi 



V47f r 



r'+^ 6 pi (r) dr . 



Using (82) in (83), one derives 



Ql = ^Ri+'Eiiuj) , 
where the following function is introduced: 

oc 

1 



(83) 



(84) 



r'"*"^ 6pi{r) dr 



I A+2 '_ e2l+l 



(85) 



and crp^(w) = crp\ j = 1,2. 



Let us rewrite the function Si (w) in the following form: 
UiM = ^^-^il-e-^), (86) 



where wq is defined as 

{i + r)-nf ^,'+^ + ie'+' 



Wo 



(87) 



Performing some transformations and accounting for 
the expression for the determinant A, one obtains the fol- 
lowing expression: 



^ , . <Pi wli - Wo 



W21 - Wo 



1-^ 



1+2 



Att \ W — Wll w — W21 ) Wll — W21 

Finally, taking into account Eqs. (32) and (39), the func- 
tion Si{w) can be represented in the following form which 
clearly demonstrates the presence of two separate modes 
characterized by the resonant frequencies uiu and uj2i- 



471" Vw^ 



Sii 



S21 



(89) 



(81) where 



S21 



3N 1 - ^'+^ top - Wll 

R^ 1-^3 y;^; - Wll 

3N 1 - ^'+^ W21 - wq 

i?| 1-^3 W2I - Wll 



(90) 
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and N is the number of delocalized electrons in the system. 

Note that for the general expression of function Ei (w) 
one should account for the damping of plasmon oscilla- 
tions and introduce the finite widths, Fu and r2i, of the 
plasmon resonances according to Eq. (62). 

Due to interaction with the external electromagnetic 
field only dipole excitations may arise in the system. There- 
fore, the case of particular interest is Z = 1. 

In the equations presented below, let us omit the sec- 
ond subscript allocated for I in the notations for plasmon 
frequencies, wi and ui2, and their widths, A and Below 
we consider only the case 1 = 1. 

The resonant frequencies of the dipole surface plas- 
mons are 



N 



N 3-p 

N 3+p _ TV 4 



(91) 



where the notation p = ^/l + 8^ is introduced. 

Calculating the values wq, w\i and W21 for Z = 1 and 

accounting for the plasmon widths, one obtains the ex- 
pression for the dipole term of the function Si{uj): 



3N(j)[ 
87rpi?i 



p+l 



i wA 



+ 



p-1 



■ i wA 



(92) 



The cross section of photoionization by a single photon 
is given by the general expression: 



0-1(0;) = ^^lmD{u)) 



(93) 



where E is the strength of the external electric field, and 
the induced dipole moment D{uj) equals to Qi (84): 



/ 4-Tr 



(94) 



Choosing the potential <pi (r) of the linearly polarized elec- 
tromagnetic wave (in the dipole approximation) in the 
form 

Mr) = ' 
and using Eq. (92), one derives 
.N 




(95) 



D{lo) = -E 



2p 



p+l 



p-1 



^uj'^ — Lol+i w A - w| -h i w A 

The imaginary part of D{lo) reads: 

N 
2^ 

(p + l) A , (p-1) A \ 



(96) 



lmD{uj) = Euj ■ 



(97) 



Therefore, one evaluates the final expression for the pho- 
toionization cross section: 



iVi A 



+ 



A^2 A \ 



where 



Ni = N 



p+l 
2p 



N2 = N 



P-1 

2p 



(98) 



(99) 



This expression clearly shows that the photoionization 
cross section is defined by the two surface plasmons. These 
plasmons were observed at the experiment [17] in inter- 
pretation made in Ref. [18]. The volume plasmon can be 
formed only when the system is exposed to a non-uniform 
external field, e.g. in collisions with charged particles. 
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